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1 . SCOPE 


1.1 GENERAL 

This specification establishes the final design of 
modifications to the ISOCLS program within the EOD-LARSYS 
system, as specified in the Requirements Document REF; 642 -n04. 





2 . 


APPLICABIE DOCUMENTS 


The following documents, of exact issue shown, form a part 
of this specification to the extent specified herein. 

e Action Document: 63-1327-169S-06 

e Requirements Specification: REF: 642-2104 

e TIRF 76-0086 
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3. SYSTEM DESCRIPTION 
3.1 HARDWARE DESCRIPTION 

N/A 


3.2 S OFTWARE DESCRIPTION 

The purpose of this modification is to: 

a. Develop a ranking capability for cluster splitting and 
cluster merging. In cluster splitting, the standard 
deviations will be ranked in descending order to insure 
that the clusters having the largest of standard devia- 
tions greater than STOMAX will be given the highest 
priority for splitting. In cluster merging, the inter- 
cluster distances will be ranked in ascending order to 
insure that the clusters having the smallest of inter- 
cluster distances less than DLMIN will be given the 
highest priority for combining. 

b. Compute and print out an Error Criterion E after each 
iteration. 

c. Develop a cluster deletion procedure after the final 

iteration. The procedure is as follows: After the final 

or statistics pass P^, there will be a test on the popu- 
lation of each cluster. If the population of any cluster 
falls below a user specified level, PMIN, plus the number 
of channels, the cluster(s) will be eliminated. There 
will be second pass ?2 reassign the pixels of the 
eliminated cluster (s) to the cluster centers of Pj^. The 
statistics pass will be there are no clusters 

deleted after Pj^, the statistics pass will be P^^. 




I 


d. Simplify the coding logical in ISODAT. To accoiitplish 
this requirement, only one procedure for clustering will 
be developed and maintained. The simulated SKIPS procedure 
will be maintained, and the coding for the standard LARSVS 
procedure will be extracted from ISODAT and SBTUP7. The 
control card OPTION SIMBRP will no longer be a valid 
control card. 


e. Expand the input capabilities of the clustering processor 
to optionally access mean values from a saved statistic 
file (cards, drum, or tape file). These mean values will 
serve as starting vectors for the cluster processor. 

The following control card inputs are required for this 
modification: 

• Cluster Ranking During Splitting and Combining 
No new control card input. 


e Error Computation 

Control card - OPTION ERCOMP 
ERCOMP is a key tliat triggers the computation of 


E «» 


1 

M 


LNCAT 

2 DN(K) 
K“1 


NOFEAT 


j-i 


* * 


oj (K) 


where DN(K) is the number of pixels in cluster K, 

LNCAT is the number of clusters, 

. NOFEAT is the number of channels used in clustering, 
Oj(K) is the variance of the jib channel of the 
' kth cluster, 

M is the total number of pixels being clustered, 
and prints the value of E after each iteration. 




• Cluster Deletion Procedure 

Control card - PMIN N 

PMIN is a user threshold value for cluster deletion. 

If the population of a cluster is less than PMIN plus 
the number of channels, the cluster is deleted after 
the final pass, P^. 

• Access Mean Values from a Saved Statistics File 

Control Cards - STATFIhE INPUT/UNIT«N,FUE-M 

STATFILE OUTPUT/UNIT-N, FILE-M 
CHANNELS STAT«Ni , . , . , DATA“Mj , . . . ,M^^ 

MODULE 

SUBCLASSES Nj,...,Nj^ 
with M and N being integer numbers. 

STATFILE INPUT/UNIT-N.FILE“M specifies the unit the 
input STAT file (starting vectors) is mounted on, and 
the number of the file to process. Defaults to self- 
starting vectors. If FILE is omitted, the default is 
File^l. 

STATFILE OUTPUT/UNIT-N, FILE«M* specifies the unit the 
output STAT file is mounted on, and the number of the 
file to create. Default values are unit«l and filesl, 


CHANNELS STAT»Nj N^ , DATA-M ^ , . . . M^ . Nj,N2,...N^ are 

channel numbers referring to the input STAT file. 

Mi,M 2 , . . .Mjjj are channel numbers referring to the MSS 
data^ape. If STATFILE INPUT/UNIT-N, FILE-M control 
card is input, STAT set of channels must be defined. 





MODUI^E triggets the Inptit of the input Module StA'i 
*cerd deck. Default«No statistics card deck* 


SUBCLASSES .. .,Njj are the subclass numbers referring 
to the input STAT file. ^ 




ISOCLS 


3.2.1 
3.2.X.1 Linkages 

ISOCLS is called from MONTOR. ISOCLS calls LARSAA routines , 
SETUP?, RDDATA, RDFILB, ISODAT, CHAIN, PRINT, PSTAPB, COVARl, 
TWRITE, and PCHSTl, UNIVAC system routines CMBRR and RWRITE 
are also utilised. 

ir 

3. 2.1.2 Interfaces 

N/A 

3. 2.1. 3 Inputs 

Input to the clustering processor ISOCLS consist of control 
cards, field/class definitions and imagery data tape. See 
SETUP? description for control card input and RDDATA for 
imagery data input. 

3. 2.1.4 Outputs 

The outputs of the ISOCLS processor are line printer cluster- 
ing results, cluster map tape and a statistics file. 

3 . 2 . 1 . 5 Storage Requirements 

Stprage used; Code-1017 (8) , Data-26475(8) . 

1*2. 1.6 Description 

ISOCLS performs a modified version of the clustering algorithm 
(ISODAT) . The program expects multispcctral scanner data in 
either the LARSYS 2 or Universal format. 

3, 2. 1.7 Detailed Flowchart 

N/A 




3.2.1. 8 Uistiag 

See Appendix A £cr program listing. 

3.2.2 SETUP7 

3. 2. 2.1 linkages 

SETUP? is called from ISOCLS. SETUP? calls NXTCHR, NUMBER, 
FLTNUM, FTND and RDMEAN. 

3. 2. 2. 2 Interfaces 

N/A 

I 

3. 2. 2. 3 Inputs 

The inputs to the SETUP? subroutine which affect this speci- 
fication are the control cards which include the keywords 
PMIN, STATFILE, CHANNELS. MODULE, SUBCLASSES and OPTION ERCOMP. 

3. 2. 2. 4 Outputs 

The output of the SETUP? subroutine is tlie input summary which 
includes cluster summary, cluster map, DAS output tape with 
cluster means (defaults) or cluster numbers. 

3.2. 2.5 Storage Requirements 

Storage used: Code“1222(8) , Data-S47(8). 

3. 2. 2. 6 Description 

SETUP? reads and analyzes all parameter values and sets de- 
faults for all control card input for the ISOCLS processor. 


N/A 




3.2.2.? Detailed Flowcharts 


1 


t 


3 . 2 . 2 , 8 Listing 

See Appendix B for program listing. 

3.2.3 ISODAT 

3. 2.3.1 Linkages 

ISODAT is called from ISOCLS. ISODAT calls UNIVAC system 
routines RROAD, RWRITE and CLOCK. LARSAA routines CLDIST 
and PRINT are also called. 

A 

3. 2. 3. 2 Interfaces 

N/A 


3 . 2 . 3 . 3 Inputs 

& 

The inputs to the ISODAT subroutine which affect this modifi- 
cation are parameter values of keywords PMIN, STATFILE, CHANNELS, 
SUBCLASSES, MODULE and OPTION ERCOMP. 

3. 2.5.4 Outputs 

The outputs of the ISODAT subroutine are the line printer 
cluster summaries at each iteration. 

3 . 2 . 3 . 5 Storage Requirements 

Storage used: Code-2010(8) , Data-614(8). 

3. 2. 3.6 Description 

All calculations of the clustering algorithms are done by 
ISODAT. 


N/A 



i 

1 

i 


i 
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Detailed Flowcharts 




3.2.4 GETSTA 

3. 2.4.1 Linkages 


GETSTA is called from ISOCLS. GETSTA calls LARSAA routine 
EXIT and Univac routines FSBSFL and CMHRR. , 

3. 2. 4. 2 Interfaces 

The module STAT file is the only interface between ISOCLS 
and other LARSAA processors. The file may be obtained from 
STAT or a previous execution of ISOCLS. 

3.2.4. 3 Inputs 

The new parameters that were added to the subroutine GETSTA 
are SUBVEC • Subclasses from input statistics file for in^cial 
means, and N0SUB2 - Number of initial means. 

3* 2.4.4 Outputs 

N/A 


3 . 2 . 4 . 5 Storage Requirements 

Storage used: Code-433(8), Data-1130 (8) . 

3 . 2 . 4 . 6 Description 

GETSTA extracts means from input MODULE STAT deck/file. The 
user may specify a subset of the channels and/or subclasses. 

3*2.4. 7 Detailed Flowcharts 

N/A 

3. 2. 4. 8 Listing 
♦ . 

See Appendix C for program listing. 




5,2.5 D£SC£N 

3.2.S.1 U.tt3Cftges 

0ESC£N is called from ISODAT. 


5. 2. 5.2 Interfacei 

N/A 

5. 2. 5. 3 litPWfeS 

The argomeats for the subroutine DESCEN are: ’ 


1 

j 


A 


SCN - An output array dimensioned by LNCAT that contains 
values in descending order. 

LNCAT - The current number of clusters. 

IPTT - An array dimensioned by LNCAT. IPTT is used as a 
pointer to the rciordered SCN array. 


3. 2. 5.4 


Output s 


N/A 


3 . 2 . 5 . 5 Storage Requirements 

N/A 

3 . 2 . 5 . 6 Description 

DESCEN ranks SCN values in descending order and IPTT is used 
as a pointer to these values. 

3. 2.5. 7 Detailed Flowcharts 

N/A 

3 . 2 . 5 . 8 Listing 

See Appendix D for program listing. 




(( 


3.2.6 PSPLIT 

3. 2.6.1 Linkages 

% 

PSPLIT is called from ISODAT. 

» 

3. 2.6. 2 Interfaces 
N/A 


3.2.6. 3 Inputs 


The arguRicnts for the subroutine PSPLIT are; 


MEANS - An input array containing the means of the last 
NMIN iteration. 


STDEV 

N 

CLD 


An input/output array containing the standaid d-- vi- 
ations. 

An input/output array containing the number of 
points in each cluster. 

An input/output array containing the intercluster 
distances. 


C - An input array for storing the imagery data, 

IPLACE - An input array for storing the cluster number each 
’ corresponding pixels belongs to. 

AVP - A temporary array for computing the means. 

MEN - An output array containing the means of the PMIN 
iteration. 


3 . 2 . 6 . 4 Outputs 
N/A 




3.2.6.5 Storage Reqmrements 
Requires S77g words of core. 

3. 2.6.6 Description 

If all clusters have a population greater than the parameter 
PMIN, PSPLIT returns to the calling routine without any re- 
assignment of the data; otherwise » all clusters with less 
than PMIN points are deleted, and the data reassigned to 
the remaining clusters. 

3. 2. 6. 7 Detailed Flowchart 

See following page. I 

I 


3. 2. 6. 8 Listing 

See Appendix E for program listing. 


I 














3.2.7 ISODAT 


3.2. 7.1 Linkages 

ISODAT is called by ISOCLS and calls the routines CLDIST, 
PRINT, PSPLIT, and DESCEN. 


3. 2. 7. 2 Interfaces 
N/A 

3. 2. 7. 3 Inputs 
N/A 

3.2. 7.4 Outputs 
N/A 


3.2. 7. 5 Storage Requirements 
Requires 30476g words of core. 

3 . 2 . 7 . 6 Description 

ISODAT perforins the clustering. The routine was modified to 
rank the standard deviations in descending order to insure 
that the clusters having the largest of standai'd deviations 
would be given the highest priority for splitting. The inter- 
cluster distances are ranked in ascending order to insure that 
the clusters having the smallest of intercluster distances 
would be given highest priority for combining, 

3.2. 7.7 Detailed Flowchart 
N/A 


3. 2. 7. .8 Listing 

See Appendix F for program listing. 








4. OPERATION 

% 

The EOD’LARSYS system is operational on the Univac 1108 
EXEC2. To exercise the capabilities implemented as describ^ed 
in this document, see the User’s Documentation BQD-yVRfYS » 
and list of added control cards in section 3.2. 


S. TEST PROCEDURE 


5.1 GENERAL 

■ " ■■ ■ ' " ■' T 

Using a 16 channel tape for checkout, vertification that th^e' 
modifications to EOD-LARSYS system were correct was substantiated 
in the following manner: 

a. The output of the cluster summary showed that the 
clusters having the largest of standard deviations 
greater than STDMAX were split and that the clusters 
having the smallest of intercluster distances less 
than DLMIN .were combined. 

b. The error criterion E was printed after each iteration. 

c. The output of the final cluster summary showed that 
after the statistics pass, there were no clusters with 
a population less than the value of PMIN plus the num- 
ber of channels. 

d. The means from the statistics tape created by the ISOCLS 
processor were used as starting vectors for clustering. 
This vertified the expansion of the input capabilities 

of the clustering processor to optionally access mean 
values from a saved statistics file. 

5.2 COMPUTER RUN 

To further vertify that the modifications were implemented 
as outlined by the Requirements Specification, Rnf; 642-2104, 
the same set of data was clustered on both ERIPS and LARSYS. 

The results are shown on the following pages. 





A 


LARSYS RUN 
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6. TEST APPROVAL SHEET 


UNIVAC 1108 LARSYS/ISOCLS MODIFICATIONS 

FOR 

LACIE PHASE III 
Job Order 71>695 
(TIRF 76-0086) 


VerificAtion Test Specification 
APPROVAL SHEET 


TEST CONDUCTOR: 

COGNIZANT ENGINEER: 

REQUIREMENTS: • 

USER: <T kT,* 




NASA TECHNICAL MONITOR: ' h\ r /?. 

MAINTENANCE 6 OPERATIONS; 

LACIE QUALITY ASSURANCE; 

QUALm' ASSURANCE! 


>URA^E: 

- -11 


REMARKS: 


DATE: • 

DATE: ^/jn/? 7 

DATE: ’;i/ich7 
DATE: ofiohy 

f f t 

DATE: /o 

DATE; 

DATE: 

DATE: f / 






6, TEST APPROVAL SHEET 


UNIVAC U08 LARSYS/ISOCLS MODIFICATIONS 

FOR 

LACIE PHASE III 
Job Order 71-695 


(TIRF 76-0086) 

Verification Test Specification 
APPROVAL SHEET 


TEST CONDUCTOR: 


COGNIZANT ENGINEER 
REQUIREMENTS: 

USER: 


^EER: 





MAINTENANCE 6 OPERATIONS: 
LACIE QUALITY ASSURANCE: 


QUALITY ASSURANCE 


REMARKS: 









NASA TECHNICAL MONITOR: OiA r J 


. ■Ti' ' 


DATE: ■ • 

DATE: A/jo /77 
DATE: -d/ 10/77 
DATE: Sfio/yj 

DATE: iQ_ 2 il 42 ? 

DATE: 

DATE: 

DATE: f f 
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